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Decoherence-induced leakage errors can couple a physical or encoded qubit to other levels, thus 
potentially damaging the qubit. They can therefore be very detrimental in quantum computation 
and require special attention. Here we present a general method for removing such errors by using 
simple decoupling and recoupling pulse sequences. The proposed gates are experimentally accessible 
in a variety of promising quantum computing proposals. 



PACS numbers: 03.67.-a,03.67.Lx,03.65.Yz 

The unit of quantum information is the qubit: an ideal- 
ized two-level system consisting of a pair of orthonormal 
quantum states. However, this idealization neglects other 
levels which are typically present and can mix with those 
defining the qubit. Such mixing, the prevention of which 
is the subject of this work, is known as "leakage" . Leak- 
age may be the result of the application of logical opera- 
tions, or induced by system-bath coupling. In the former 
case, a rather general solution was proposed in Q. Here 
we are interested in decoherence-induced leakage. E.g., 
in the ion-trap QC proposal the two-level approximation 
may break down and spontaneous transitions may leak 
population out of those levels that represent the qubit 
in the ion [0. This is part of a more general problem: 
quantum computation (QC) depends on reliable compo- 
nents and a high degree of isolation from a noisy environ- 
ment. When these conditions are satisfied, it is known 
that it is possible to stabilize a quantum computer us- 
ing an encoding of a "logical qubit" into several phys- 
ical qubits. Methods which profitably exploit such an 
encoding are, e.g., (closed-loop) quantum error correct- 
ing codes (QECC) || || and (open-loop) decoherence- 
free subspaces or subsystems (DFS) f| 0, |J. The 
logical qubits of these codes can also undergo leakage 
errors, which are particularly serious: by mixing states 
from within the code and outside the code space, leakage 
completely invalidates the encoding. A simple procedure 
to detect and correct leakage, which can be incorporated 
into a fault-tolerant QECC circuit, was given in ||. This 
scheme is, however, not necessarily compatible with all 
encodings || . Here we present a universal, open-loop so- 
lution to leakage elimination, which makes use of fast and 
strong "bang-bang" (BB) pulses [[[o[ [n]. We first give 
a general scheme for protecting qubits (whether encoded 
or physical) from leakage errors using an efficient pulse 
sequence. Then we illustrate the general result with ex- 
amples taken from a variety of promising QC proposals. 
Particularly important is the fact that our scheme is ex- 
perimentally feasible in these examples, in the sense that 
we only make use of the naturally available interactions. 

Universal leakage elimination operator. — Here we give 
a general, existential argument for eliminating all leakage 



errors on encoded or physical qubits. Suppose that n two- 
level systems (e.g., electron spins in quantum dots | fl2| ) 
are used to encode one logical qubit, or that a TV-level 
Hilbert space Hjv supports a two-dimensional physical 
qubit subspace (e.g., hyperfine energy levels of an ion 
||[). Let us arrange the basis vectors {\n)}^=o of Hk so 
that |0) and jl) represent the (physical or encoded) qubit 
states (N = 2™ for the encoded case). We refer to this 
as the "ordered basis" . In this basis we can classify all 
system operators as follows: 




where B and C are 2 x 2 and (N - 2) x (N - 2) blocks 
respectively, and £>, F are 2 x (TV - 2) , (TV - 2) x 2 blocks. 
Operators of type E represent logical operations, i.e., 
they act entirely within the qubit subspace. £ ,J - oper- 
ators, on the other hand, have no effect on the qubit as 
they act entirely outside the qubit subspace. Finally, L 
represents the leakage operators. The total system-bath 
Hamiltonian can be written as Hsb = He + H E ± + Hr, , 
where He {H e ± , Hr,) is a linear combination of elements 
of the set E (E 1 -, L), tensored with bath operators. Now 
consider 

where the blocks have the same dimensions as in Eq. ([!]) 
and 4> is an overall phase. This operator satisfies 
{R L ,L} = 0, while [R L ,E] = [RliE 1 -] = 0. Using a BB 
parity-kick sequence |l0[ it follows that Rr, is a leakage- 
elimination operator (LEO): 

lim ( e - lHsBt/n RU- lHsBt/n R L ) n = e -iH E t e - t H E± t 

n — >-oo 

(3) 

In practice one takes n = 1 and makes the time t very 
small compared to the bath corr elation time @. Eq. (§) 
then holds to order t 2 , and implies that one intersperses 
periods of free evolution for time t with Rl,R\ pulses 
which are so strong that Hsb is negligible during these 
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BB pulses. The term e~ lH E ±t in Eq. (||) has no effect on 
the qubit subspace. The term e~ lHEt may result in logi- 
cal errors, which will have to be treated by other meth- 
ods, e.g., concatenation with a QECC || [|, [H| or ad- 
ditional BB pulses jll|, [l5| . Note that since Rl commutes 
with the logical operations, they can be performed at the 
same time, i.e., our leakage elimination procedure is fully 
compatible with universal QC. We now give a procedure 
for generating LEOs from a controllable system Hamil- 
tonian H$ acting for a time r, i.e., Rl = exp(— iHsr). 
From Eq. (|^) it follows that H$ must act as a projec- 
tion operator P onto the qubit subspace. Furthermore, 
r must be chosen so that Rl acts as —I in the qubit 
subspace. A general choice is 

= exp (±iTrn ■ aP) (4) 

where a denotes the vector of Pauli matrices, which we 
refer to as logical X, Y, Z operations, and n is a real unit 
vector. This is a valid LEO since exp {—mri ■ a) expresses 
a 2-7T rotation about the axis n on the qubit Bloch sphere, 
upon which the qubit state acquires a minus sign. A use- 
ful example is r = ir and Hs = |0)(0| + |1)(1|, which 
is a projector onto the qubit subspace and acts as iden- 
tity there. This example generalizes immediately to d- 
dimensional qudits: 

d-i 

i?^ d(<i) =exp(± l ^^|fc}(fc|). (5) 

k=0 

Now let us consider leakage prevention on a code sub- 
space of K logical qubits, each supported by n physical 
qubits. In analogy to R L we can construct a general 
LEO as follows. Let Si be a single-qubit logical (unitary) 
operation on the i th (encoded or physical) qubit, and let 
Pi be a projection on the code subspace of that qubit. 
Then 

Rl [K) = exp(±i7r ®f =1 S t Pi) (6) 

is a valid LEO. 

Proof: We can always rotate the Bloch sphere of a 
qubit so that each Si is independently transformed into 
Zf. Si — UiZillJ (where f/j is an appropriate single- 
qubit unitary); ®fL 1 Zi is a diagonal matrix of ±1, so 
exp(-i7r <Z>f =1 Zi) = -I. Thus R E L (K) = exp(-i7r 

UiZiUjPi) = (®? =1 Ui)exp(-™®? =1 ZiPi){®?=iUt) = 
— I on the code subspace, and (because of the Pj) = I on 
the orthogonal complement. QED. 

Sometimes we shall be able to construct single-qubit 
logical operators which are automatically projectors on 
that qubit's subspace. We refer to such operators as 
"canonical" . We are now ready to apply these consid- 
erations to a number of promising QC proposals. 

Example 1. — As a simple first example, consider phys- 
ical qubits (without encoding), such as electrons on liq- 
uid helium [116), or an electron-spin qubit in quantum 



dots |12j, h7|, or a nuclear-spin qubit in donor atoms 
in silicon |18|. In those potential well at each 

site traps one fermion. Usually, the ground and first 
excited state are taken as a qubit for a given site: 
\k) = ctjvac), where ct is a fermionic creation operator 
for level k — 0,1. Let n& = cj.Cfc be the fermion num- 
ber operator. The logical operations for this qubit are 
E = {X = cjci + cjci, Y = i(c\c - cjci), Z = n Q - n x } 
whose elements satisfy su(2) commutation relations. In 
this case, a general linear Hamiltonian which includes 
hopping terms, H S b = Y,k7=o a kic\ci, where au includes 
parameters and bath operators, and k, I denote all elec- 
tron states, can leak the qubit states k = 0, 1 into any 
of the other states. Using parity-kicks, we can elimi- 
nate this leakage in terms of the LEO [recall Eq. (|J)]: 
R 1 ^ 1 ^ = exp[±i7r(no + Ui)]. This LEO is implemented 
simply by controlling on-site energies. 

Let us now generalize this to K qubits. The states 
of the i th qubit are \k)i = ct(i)|vac). States outside of 
the code subspace contain at least one creation operator 
c\(i) with k > 2. A logical Z operator on the i th qubit 
is Zi = rio(«) — ni(i), which is canonical. It follows from 
Eq. (§) that an LEO is: 

(i?i (K) )fcrm = exp(±t7T(Z 1 Z 2 ■ ■ ■ Z K )). (7) 

The term ®^ =1 Zi involves a many-body interaction which 
is not naturally available. However, it can be constructed 
from available interactions as follows: Let us assume that 
the interaction between neighboring sites i, j contains a 
controllable ZiZj term (in reality such control may have 
to be obtained indirectly, e.g., by controlling an XiXj + 
YiYj term, as shown in |l9[, and as discussed in more 
detail below). We note the following useful "conjugation 
by tt/4" formula 0: 

rp i8B _ if A iBB —i%A iBUAB) r Q \ 

±A°e =e 4 e e 4 =e y ', (8) 

which holds provided {A, B} — and A 2 = B 2 = I. 
Using this we can efficiently generate long-range interac- 
tions by alternately switching interactions A, B on/off. 
E.g., 

Ty 2 o [T Z2Z3 o (T X2 o e^)] = e «™> z >, (9) 

which can, in turn, be used to generate e * ez i- z 2Z 3 ^4 ) e £ a 
Using this recursive process, the implementation of the 
LEO {R^ K ^ ) term takes O (K) steps. Fig. [I] shows a circuit 
for the 4-qubit case. 

Finally, we note that we can also treat bosonic sys- 
tems, such as the linear optical QC proposal ^0|. In 
this case, a qubit is encoded into two modes. The first 
qubit has states |0)i = 6||vac) and |l)i = foj I vac) , and 
the second qubit is |0)2 = 6g|vac) and 1 1 ) 2 = ^||vac), 
where b\ are bosonic creation operators. Encoded two- 
qubit states are |00) = b\b\\v&c), |01) = b\b\\v&c), 1 10) = 
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FIG. 1: Quantum circuit of the LEO for 4 physical (encoded) 
qubits in Example 1 (3). Ua = e z ~z , where A stands for 
X i% Yi,Zi,ZiZi (XuY^ZuZiZj) for Example 1 (3). 

6^3 1 vac) and |11) = b^lvac). But the linear optical 
Hamiltonian H = l=1 akib\b\, contains beam-splitter 
terms like b\b^ and 6^3 > which can cause leakage into 
states such as &|&2l vac ) or b^\va,c). By using the LEO 
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id(i) 



exp 



we can eliminate the 



±m{b\b l + b\b 2 ) 

leakage terms. This LEO can be implemented simply us- 
ing a linear optical phaseshifter. However, generalizing 
this LEO to multiple encoded qubits requires a photon- 
photon interaction, which is not readily available. 

Example 3. — A substantial number of promising solid- 
state QC proposals, e.g. JjJ, [H], |J , are governed 
by effective isotropic and anisotropic exchange interac- 
tions, which quite generally, can be written as 



He* — ^ JfjXiXj + JfjYiYj + J^Z l Z j 



(10) 



The encoding |0) L = |01), |1) L = 1 10) (using two physical 
qubits per logical qubit) is highly compatible with H ex , 
in the sense that universal QC can be performed by con- 
trolling the single parameter J?j in the Heisenberg (J?- = 
4 = J?j), XXZ (J* = ±4 ^ J?.), and XY (Jg = % 
J?j = 0) instances of i? C x, provided there is a Zeeman 
splitting that distinguishes single-qubit Zi terms. This 
is done using the "encoded selective recoupling" method 
|L9||. Furthermore, the {|01), |10)} encoding is a DFS for 
collective dephasing (where the bath couples only to sys- 
tem Z 2 i-i + Z 2i operators) (|, [l|, ^2|. A set of logical 
operations on the code is E = {X\ = {X\X 2 + YiY 2 )/2, 
Y x = {X 2 Y 1 -Y 2 X 1 )/2,Z 1 = (Zi-Z 2 )/2}. OnlytheAr 
term is assumed to be directly controllable (by manipula- 
tion of Jf 2 ), while the Z\ term can be turned on/off using 
recoupling |19[ . The Y\ term can then be reached in a 
few steps: e~ l0Yl — T x - i o e~ l9Zl . The leakage errors are 
due to system-bath interactions where the system terms 
include any of Xi,Yj, Xi Zj and Yi Zj , since as is eas- 
ily seen, such terms do not preserve the {|01), |10)} code 
subspace. As pointed out first in Jig ], the LEO can be 
expressed as \R L ) 2 -dfs = exp(i7rAi) = ZiZ 2 , which 
means that it is implementable using just the controllable 
Jf 2 parameter in the instances of H ex mentioned above. 
This form for (-R^^-dfs is an instance of Eq. (H), 



with n — x. Note that, in agreement with our general 
comments above, (R L ) 2 -dfs commutes with every el- 
ement of B 2 _dfSj meaning that logical operations can 
be performed on the encoded subspace while eliminating 
leakage. 

Next we now show how to efficiently eliminate leakage 
in this case on an arbitrary number of encoded qubits. 
The m th logical qubit is encoded as |0l) to = |0 2m _il 2 m), 
li)m = |l2m-i0 2m ). The logical Z operator is Z m = 
(Z 2m -\ — Z 2m )/2, and is canonical. It follows from 
Eq. (|) that a valid LEO is: 



{-n-L J2-DFS 



exp(±i7rZiZ 2 ...Z k)- 



(11) 



The next question is how to efficiently construct such an 
operator. The term Z\Z 2 = j(ZiZ^ + Z 2 Z^ — Z 2 Z?, — 
Z\Z±) contains next and second- next nearest-neighbor 
interactions. Using "conjugation by tt/4", they can all 
be generated using only nearest-neighbor interactions in 
terms of the relation 



o i0 ZiZi-^-2 rp 



= T Xi+lXi+2+Yi+lYi+2 oe iez * z ^ (12) 

where, in accordance with Jl9| , we have only assumed 
controllability of the XY interaction term Xi + iXi +2 + 
Yi+{Yi-\-2. At this point we can use the recursive con- 
struction of Eq. (g) again, by replacing X, Y, Z there 
by their encoded counterparts. Doing so takes ZiZi + \ 
to ZiZi + \Zi +2 , etc., and will again efficiently construct 

the LEO (R^} K ^) 2 -dfs, i-e., using O(K) steps. An ex- 
ample of this for 4 encoded qubits is shown in Fig. [l]. It 
is interesting to contrast the linear scaling of this leak- 
age elimination procedure with general error elimination 
using BB pulses. As shown in [n], without additional 
symmetry assumptions restricting the order of coupling 
terms in the Hamiltonian Eq. (||), the BB procedure, if 
used to eliminate all errors, requires a number of pulses 
that is exponential in N. 

Example 4- — Collective decoherence is a system-bath 
interaction that obeys full qubit permutation symmetry: 
tfColLDec = J2 a=x y z (£ i0f ) ® B ° , where a? are the 
Pauli matrices and B a are bath operators || Q, [t) - This 
situation can be created from an arbitrary linear system- 
bath coupling = J2i Bi, where a l = (X t , Yi , Zi) 
and Bi are bath operators, using a BB symmetriza- 
tion pulse-sequence that employs only the Heisenberg ex- 
change interaction J23|. The shortest DFS (or "noiseless 
subsystem" ) encoding that protects a single logical qubit 
against collective decoherence uses 3 physical qubits M . 
In |2^] it was shown that one can perform universal QC 
on this DFS, again using only the Heisenberg interaction. 
To explain the encoding, note that the Hilbert space of 3 
spin 1/2 's has total-spin S = \{a\ + o 2 + 173) and splits 
into two 5* = 1/2 subspaces (denoted A = 0, 1), and aS = 
3/2 subspace. The states can be labeled as \S, A, <!5 2 ), and 
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a|l/2, 0, 1/2} + /3|l/2, 0, —1/2), which again is a procedure that scales as O(K). 



the DFS qubit is \0 L 
\1 L ) = a\l/2, 1, 1/2} + p\l/2, 1, -1/2}, \a\ 2 + \p\ 2 = 1, 
i.e., the encoding is into the degeneracy of the two 
5 = 1/2 subspaccs [|[ 22 . Collective errors can change 
the a, (3 coefficients, but have the same effect on the 
\0l), |1l} states, which is why this encoding is a DFS. 
If, however, we also consider bilinear system-bath cou- 
pling Hfl 



^2i<j Ea,)3={i,ii,2} 9if (J ?' J 



B°f, where 

is a rank-2 tensor, then the symmetrization proce- 
dure of [^3), that prepares collective decoherence condi- 
tions, will not work. In this case we must consider the 
possibility of leakage. The bilinear term gtf' J t <J ^ can be 
decomposed into (i) a scalar gai ■ <jj , which has the effect 
of logical errors E; (ii) two rank-1 tensors (3- (Si x <jj) 
and (di ■ 7) (ffj ■ 7), which can couple between 5=1/2 
states, and can couple them to 5 = 3/2 states. Note that 
this also applies to imperfect symmetrization at the level 
of a linear system-bath Hamiltonian Hg£ . Thus we see 
that the 5 = 3/2 subspace acts as a source for leakage, 
and that there is also the possibility of (non-collective) 
errors [both from (ii)] which do not have the same effect 
on the |0l), |1l) states. We defer an analysis of the latter 
"5 = 1/2 — > 1/2" errors to a separate publication pi] , 
but we note that they can be suppressed using techniques 
similar to those we discuss next. 

An open-loop leakage correction circuit for this DFS, 
that once more uses only the Heisenberg interaction, was 
given m @. There the DFS qubit was defined to be 
|0 L } = |l/2, 0,1/2), |U) = |l/2,l,l/2} and transitions 
to any of the other 6 states were considered as leakage 
(this includes errors caused by collective decoherence, 
which arc normally avoided by a DFS encoding). Here 
we add another element to this picture of the Heisenberg 
interaction as an enabler of universal, fault-tolerant QC, 
by showing that it can also provide an LEO. The impor- 
tance of Heisenberg-only QC, as pointed out in [p5| , is 
in the relative ease of manipulating this interaction in a 
number of the most promising solid-state QC proposals 
IfL^liql- Now, as shown in [|||, X — — -1x3 -a 2 -a 3 ) 
acts as a logical X on the DFS qubit defined above, 
and annihilates the 5 = 3/2 states (i.e., it is canonical). 
Therefore, using Eq. (||), (i?^ 1 ^3-DFS = exp(±z7rX) is 
a Heisenberg-only LEO for a single 3-qubit DFS which 
eliminates transitions to the 5 = 3/2 subspace. An LEO 
for the if-qubit case is then, from Eq. (H): 



{Rf K) ) 



-dfs = exp(±i7r(XiX 2 • ■ • X K )). 



(13) 



To generate this LEO from available interactions we 
use a procedure similar to Eq. (||). First, note that 
TV oTy o e l9Xl = e* 9 * 1 ^ 2 . Efficient schemes 

for generating ZiZj were given in [ p5[ , while Zi,Yi are 
directly obtainable from the Heisenberg interaction [B2|. 



The recursive construction of (R L 
ceeds using Ty. ° o ° e 



E(iC), 



Conclusions. — Decoherence-induced leakage from the 
logical space of (physical or encoded) qubits is a severe 
source of errors for quantum computation. We have 
shown how to efficiently and universally eliminate such 
errors using sequences of "bang-bang" pulses. These 
pulses can be applied at the same time as logical op- 
erations, so that leakage elimination can be performed in 
conjunction with universal quantum computation. Ap- 
plications to a variety of promising quantum computing 
proposals were discussed, and leakage elimination meth- 
ods were presented that are directly applicable using only 
experimentally available interactions. 

The present study was sponsored by the DARPA- 
QuIST program (managed by AFOSR under agreement 
No. F49620-01-1-0468), by PRO, PREA, and the Con- 
naught Fund (to D.A.L.). We thank Dr. S. Schneider for 
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